Analysis of broadband microwave conductivity and permittivity measurements of 

semiconducting materials 
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We perform broadband phase sensitive measurements of the reflection coefficient from 45 MHz up 
to 20 GHz employing a vector network analyzer with a 2.4 mm coaxial sensor which is terminated by 
the sample under test. While the material parameters (conductivity and permittivity) can be easily 
extracted from the obtained impedance data if the sample is metallic, no direct solution is possible 
if the material under investigation is an insulator. Focusing on doped semiconductors with largely 
varying conductivity, here we present a closed calibration and evaluation procedure for frequencies 
up to 5 GHz, based on the rigorous solution for the electromagnetic field distribution inside the 
sample combined with the variational principle; basically no limiting assumptions are necessary. A 
simple static model based on the electric current distribution proves to yield the same frequency 
dependence of the complex conductivity up to 1 GHz. After a critical discussion we apply the 
developed method to the hopping transport in Si:P at temperature down to f K. 
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I. INTRODUCTION 

The rapid development of communication as well as 
industrial and medicine technologies demands accurate 
characterization of components at ever increasing fre- 
quencies, i.e. beyond the radio frequency range; this be- 
comes in particular relevant for insulating and semicon- 
ducting materials employed in electronic devices^. Beside 
the high technological relevance of doped semiconductors 
they are also subject to intense fundamental research 
concerning disordered electronic systems with electron- 
electron correlations^. 

On a macroscopic scale and under steady-state condi- 
tions, the interaction of a material with electric field is de- 
termined by its conductivity and dielectric permittivity. 
The desired broadband characterization of those para- 
meters becomes more challenging with rising frequencies 
because losses and spatial variation of current and vol- 
tage then gain importance. 

The materials characterization up to the MHz range 
turns out to be comparably simple: the voltage drop is 
measured when a current passes homogeneously through 
the specimen; lock-in technique allows for the determina- 
tion of the complex response. As soon as the GHz range is 
approached, the wavelength becomes comparable to the 
leads and specimen, waveguides have to be utilized and 
reflection or transmission coefficients are measured. In 
this spectral range a vector network analyzer is a suited 
and powerful tool. 

While standard circuit theory applies to radio frequen- 
cies, in the microwave range the wavelength becomes as 
short as a few millimeters and thus a careful treatment of 
the electromagnetic field distribution within the sample 
is necessary to obtain the material parameters from the 
impedance data gained by the measurement. Whereas 
the evaluation is straightforward for metallic samples un- 
der investigational^, only some approximate solutions 
and models have been developed in the past to treat the 
dielectric materials^ - — , considering in most cases liquids 
or soft matter at ambient conditions. In the course of 



investigating the dynamical conductivity of doped semi- 
conductors at low temperatures, we revisited the existing 
methods and elaborated a reliable evaluation procedure 
with optimized theoretical and experimental complexity. 
We are now able to obtain both the real and imaginary 
parts of the conductivity a = o\ + i<Ji (or the dielectric 
function e = e\ +162) by a broadband measurement from 
0.1 GHz to 5 GHz. (Throughout this text, e is the com- 
plex dielectric function of the material, relative to the 
free space permittivity erj.) 

As an application, the conductivity and permittivity of 
insulating Si:P at temperatures down to 1 K in a broad 
frequency and donor-concentration range have been in- 
vestigated. In particular, those measurements provide 
new information on the influence of electronic correla- 
tions on the hopping transport in Si:P as well as on the 
critical scaling of its dielectric constant at the metal- 
insulator transition driven by the doping concentration. 

We start in Sec. lU with the broadband microwave mea- 
surement of a general solid sample terminating the coa- 
xial line. Following a brief survey of the well-known eva- 
luation methods for metallic materials, the problems arc 
formulated which arise for samples with non-metallic con- 
ductivity. A simple static model is presented in Sec. IIIII 
for the current distribution inside a semiconducting sam- 
ple, valid in the low-frequency range up to 1 GHz. In 
Sections IIVI and [V] it is followed by a rigorous solution 
to the problem making no severe simplifications. Finally, 
the method is applied to the hopping transport in Si:P 
in Section IVTl 



II. BROADBAND MICROWAVE 
CONDUCTIVITY AND PERMITTIVITY 
MEASUREMENTS 

The experimental arrangement for complex permitti- 
vity measurements, where the material under investiga- 
tion is placed at the aperture of a coaxial probe, allows 
for a broadband phase sensitive measurement of the re- 
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FIG. 1: (color online) Schematics of a broadband microwave 
spectrometer that employs a vector network analyzer and co- 
axial waveguides. A solid sample terminates the otherwise 
open-ended coaxial line (after Ref. @). 




current flow 



FIG. 2: (color online) Current flow (dark magenta area) in 
a metallic sample (bright yellow) with a contact mask (light 
grey) on top of it; in our case 2a — 0.6 mm and 2b = 1.75 
mm. The relevant limiting cases a) d <C 5 and b) d S> 5 lead 
to a simple evaluation procedure (S < l[im is the skin depth). 



flection coefficient using a vector network analyzer. In 
case of solid matter some smart jig is needed to press the 
flat sample surface against the sensor^, as depicted in 
Fig. [1] Metallic contacts (gold or aluminium) are usually 
evaporated on top of the solid specimen that match the 
inner and outer conductors of the coaxial probe in size. 
These contacts provide a proper electrical connection be- 
tween the sample material and the probe, and they define 
the geometry of the sample surface exposed to the signal. 
The sensor can be inserted into a temperature-controlled 
environment and, for a given probe size, there is a useable 
frequency range as broad as two orders of magnitude. 

The problem of extracting the interesting material pa- 
rameters from the measured reflection coefficient data 
can be solved in two steps: first, one needs to obtain the 
complex sample impedance from the measured reflection 
coefficient (S parameter), and second, the complex mate- 
rial properties have to be calculated from the impedance. 



A. Evaluation procedure for metallic samples 

The evaluation of metallic samples has been developed 
and experimentally tested in the recent years^: 

The first task is to obtain the complex sample 
impedance Z from the reflection coefficient 5n iTO , mea- 
sured by the test set of the network analyzer, like the 
HP 8510. The general error model for a reflection 
measurement^ results in the following relation: 



Sn 



Ed 



ErSi 



1 — EsSn 



(1) 



between the measured S parameter and the actual re- 
flection coefficient Sn of the sample. The three inde- 
pendent complex values Er, E$ and Ed comprise the 
contribution of the microwave line. To determine those, 
measurements of three independent calibration samples 
with known actual reflection coefficients Sn as functions 
of frequency and temperature are required. We use bulk 
aluminium samples as short, teflon samples as open and 
thin metallic NiCr films as load standards. 



The sample impedance Z then follows directly viai£: 

(2) 



l-Su 

where Zq = 50 f2 is the characteristic impedance of the 
microwave line. 

The way to extract the complex electric conductivity 
a of a metallic sample from its complex impedance Z is 
straightforward, if the thickness of the specimen d either 
significantly exceeds the skin depth 5 or vice versa: 

• d <C S. (Fig. [5^) In case of a thin film evaporated on 
an insulating substrate, the electric field strength 
stays nearly constant throughout the whole film 
thickness d. The relation between the conducti- 
vity a and Z depends only on the geometry of the 
contacts. For the ring of inner radius a and outer 
radius b between the contacts, it reads: 



1 ln{6/a} 
a 2nd 



(3) 



d 3> S. (Fig. [2b) For typical microwave frequencies, 
the electromagnetic wave in a thick metallic sam- 
ple is significantly damped already at the depth 
of 1 /im by the skin effect. Hence, the interaction 
with the incident wave takes place in a thin layer at 
the sample surface. Boundary effects at the edges 
of the relatively broad contact area (cf. Fig. [5]) are 
negligible and the concept of the surface impedance 
Zs based on the assumption of plane wave propa- 
gation works well 4 . In case of a ring with inner and 
outer radii a and b, the formula to extract the con- 
ductivity from the measured impedance Z reads^: 
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with the magnetic permeability of vacuum /j,q and 
the free space permittivity Eq. 
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B. Formulation of the problem for semiconducting 
materials 

For materials with non-metallic conductivity, none of 
the simple assumptions valid for metallic samples hold 
and both steps of the evaluation procedure contain addi- 
tional challenge: 

1. For the calibration of the microwave line the open 
standard is as significant as the short and the load 
standards, when insulating samples are measured. 
It acts as a complex capacitor and the dc assump- 
tion S°f cn = 1, corresponding to Z opcn — oo, is 
not sufficient at higher frequencies. The correct 
frequency dependence of the reflection coefficient 
S°f cn is indispensable, as shown in Sec. [V] 
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FIG. 3: (color online) Sketch of a semiconducting sample with 
the aluminium contact-layer (light grey) on top. 2a = 0.6 
mm, 2b = 1.75 mm, corresponding to a 2.4 mm microwave 
adapter, c = 5 mm, d = 2 mm 



2. The electromagnetic wave penetrates deep into a 
non-metallic sample because - in contrast to a 
metal - the real part of the dielectric permittivity 
is positive. Field decay in the sample is only caused 
by the geometry and determined by the electric 
contacts (and, only to some minor degree, by the 
low absorption due to the imaginary part of e). In 
case of a 2.4 mm coaxial probe, the electric field 
strength falls below 1 % of its original value at the 
sample surface only after penetrating more than 
2.5 mm, for frequencies up to 5 GHz and relative 
dielectric constant up to 50. With other words, the 
penetration depth of the electric field is of the order 
of the contact area dimensions (compare Fig. [3]) and 
increases rapidly with rising frequency and permit- 
tivity. Thus, in contrast to the metals, the spatial 
field distribution which forms inside an insulating 
semiconductor is significantly different from that 
of a plane wave and it depends on frequency. Its 
knowledge is essential to extract the complex con- 
ductivity a (or permittivity e) of such a sample 
from its complex impedance Z . Integral equations, 
implying the accurate solution for the electromag- 
netic field in the sample, cannot be directly solved 
for e. Hence, approximations are required or sim- 
plified models need to be developed, with a limited 
range of validity for frequency and permittivity. 

In the next three Sections we treat the problem of 
extracting the material parameters of a semiconducting 
sample from the impedance data. We first suggest in 
Sec. lIIIl a simple static model which yields already a good 
approximation. In Sec. II VI we make use of the results by 
Levine and Papas^i and Misraii to relate the complex 
impedance of a semiconducting sample to its complex 
permittivity. In the following, the rigorous electromag- 
netic field distribution from Sec. |IV] is used in Sec. |V] to 
determine the frequency dependence of the reflection co- 
efficient .Sqf 611 of the open standard for the calibration 
step mentioned above. 



III. A STATIC MODEL FOR THE RELATION 
BETWEEN THE SAMPLE IMPEDANCE AND 
THE COMPLEX CONDUCTIVITY 

As a first approach to extract the conductivity of a 
low-loss semiconducting sample from its impedance, we 
developed a simple static model for the current distribu- 
tion in a sample, based on the following assumptions: 



1. The response is local: 

j(r) = a ■ E{r) , 



(5) 



where j is the electric current density and E is the 
electric field vector. This assumption implies that 
the electric field does not vary significantly at the 
distances of the mean free path £, which in the case 
of hopping transport is the mean separation of the 
hopping partners. 

2. The dependence of the electric field E(z,p) on the 
cylindrical coordinates z and p can be accounted for 
separately. There is no dependence on the angular 
coordinate tp due to the radial symmetry of the 
problem. 

3. Inside the coaxial line the principal TEM mode is 
excited; thus, only the radial component of the elec- 
tric field E(p) exists. The Gauss theorem yields 
E(p) = const./ p. With the voltage U between the 
mask contacts of radii a and b (Figs. [3] and [3J it 
follows: 
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(6) 



4. As far as the z dependence of the electric field 
strength is concerned, we assume that the field is 
concentrated at the surface and gets weaker for fur- 
ther depth because the path length for the corres- 
ponding current element dl increases. 
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FIG. 4: (color online) Geometry of the current distribution in 
a semiconducting sample of thickness d, with metallic contacts 
at distance 1(0) = (6 — a) (cf. Fig. [3]) as assumed in the static 
model. 



In order to calculate the total current flowing 
through a sample, we have chosen the cross sec- 
tion of the sample at the mid-distance between the 
Al-contacts (Fig. [3]). The single current line is ap- 
proximated by a triangle shape with the apex at 
the mid-distance cross section. Hence, we consider 
each infinitesimal current line at its lowest point 
designated by the z coordinate and assume for the 
corresponding electric field E(p, z) to be recipro- 
cally proportional to the length of the current line 
in order to keep the voltage U constant: 



E(p,z) 



E (p) ■ 77T • 
l(z) 



(7) 



where l(z) = 1(0) y/1 + [2z/l(0)f 



Now that E(p, z) is constructed, we can calculate the 
total current / flowing through a semiconducting sample 
using Eqs. ([5])- (JTJ) . The integral is taken over the entire 
mid-distance cross section (see Fig. [4]) with the infinite- 
simal element dq 
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Uo)J 



where 1(0) = b — a (Fig. 0|. Thus, we have obtained a 
relation between the complex impedance Z = U/I at the 
sample surface and the complex conductivity a: 



1 

a= z 



ln{6/a} 



n(b — a) arctan{2d/(& — a)} 



(8) 



In the limiting case of a thin conducting film there is a 
simple geometrical relation between the conductivity and 
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FIG. 5: (color online) Frequency dependent conductivity of 
Si:P samples with relative donor concentration n r — 0.56 and 
0.9 measured by a vector network analyzer at T = 1.1 K. The 
solid lines represent the conductivity obtained by the evalu- 
ation procedure described in Sec. IIVI The dash-dotted lines 
result from evaluating the complex impedance data using the 
formula iJHJ of the static model. The dashed lines stem from 
the primitive assumption S°f cn = 1 for the open calibration 
standard. Concerning the conductivity, all the three lines for 
either concentration n r are virtually identical, with only slight 
effects above 1 GHz. 



the impedance because the z dependence of the electric 
field and the boundary effects can be neglected. This 
allows us to check the above formula in the limit d/l(0) — > 
0, where we recover the Eq. (TjJ|). 

The static model has been applied to analyze frequency 
dependent impedance measurements of Si:P at low tem- 
peratures. As demonstrated by the dash-dotted lines in 
Figs. [5] and [U the results of Eq. (jSJ) agree very well with 
the rigorous solution outlined in the following Sections. 
Deviations can be noticed only above 1 GHz in the di- 
electric constant. 



IV. RELATION BETWEEN THE COMPLEX 

IMPEDANCE AND THE COMPLEX 
PERMITTIVITY OF A SEMICONDUCTING 
SAMPLE 

A. General considerations 

A large variety of methods have been developed in 
the past to extract the properties of low-loss and lossy 
dielectrics from a reflection coefficient measurement in 
the radio-frequency and microwave range. Simple inge- 
nious models and analytical solutions, that are valid in 
a limited parameter range, have been suggested; time- 
consuming but arbitrarily precise numerical approaches 
have been treated depending on the specific practical 
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FIG. 6: (color online) Typical permittivity spectra of Si:P- 
samples measured by a vector network analyzer at T = 1.1 K. 
The solid lines correspond to the evaluation procedure de- 
scribed in Sec. IIVI The dash-dotted lines result from evalu- 
ating the complex impedance data using the formula (JSj of 
the static model. The dashed lines stem from the primitive 
assumption S^f 011 = 1 for the open calibration standard. 



goals. A comprehensive list of references is available in 
review articles like Ref. [I?]- 

In studies of soft and liquid materials, the coaxial 
probe was frequently modelled as an equivalent circuit 
consisting of several fringe-field capacitors in the lumped- 
element approac h 7 ' 12 i 14 i 15 . In Ref. [12 a comprehensive, 
detailed and critical revision of this method can be found. 
The most striking point is the strong dependence of the 
model capacitances on the permittivity of the material 
that terminates the coaxial line, thus the approach is 
limited to specimens with dielectric properties close to 
those of the reference materials available. 

Here we consider a convenient analytical way to extract 
the complex conductivity a from the sample impedance 
Z based on the works of Levine and Papas^i and Misraii. 
The method is valid at least up to 5 GHz for the 2.4 mm 
probe and relative dielectric constants up to 50; an ex- 
tension to higher frequencies is possible with certain nu- 
merical procedures added. As an intermediate result, 
there is an integral expression for the sample admittance 
Y = l/Z as a. function of the material dielectric func- 
tion e, which is well suited to determine the frequency 
dependence of the open calibration standard in a closed 
manner (Sees. iHBl and W} . The theoretical expressions 
for the electromagnetic field on both sides of the sensor 
aperture^ are rewritten for the case of a medium with an 
arbitrary complex permittivity in the sample half-space. 
For the parameter range considered here^, we regard 
the variational principle applied by Levine and Papas as 
preferable to the precise but time-consuming numerical 
method of point-matching proposed by Mosig^^. 
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FIG. 7: (color online) The plane z = constitutes the inter- 
face between the coaxial wave guide (z < 0) and the sample 
space (z > 0). It serves as a reference plane, at which the re- 
flection coefficient Su, the impedance Z and the admittance 
Y of the sample are defined. 



B. Solution in form of an integral equation 

In the following, the coaxial wave guide with a center 
conductor of radius a and an outer conductor of radius 
b is terminated by an infinite-plane conducting flange at 
z = (Fig. [7|). Choosing the dimensions a and b of 
the coaxial line to be small enough, the assumption of 
a single propagating mode (the principal TEM mode) in 
the coaxial region is justified in the covered frequency 
range. This system is amenable to a detailed theoreti- 
cal analysis^ that yields the electromagnetic field distri- 
bution in the half-space z > and a relation between 
the aperture admittance Y (current-to-voltage ratio at 
2 = 0) and the complex wave vector k of the free space 
z > 0. 

We may assume an insulating sample to fill the half- 
space z > by choosing its finite dimensions to be large 
enough for the electromagnetic field strength so that the 
sample boundaries are negligible (cf. Sec. Ill Bp . The 
results for the free space thus transform into a relation 
between the admittance Y measured at the sample sur- 
face z — and its complex dielectric function sphere, a 
harmonic time dependence exp{— iujt} is assumed for the 
electromagnetic field^): 



Y 



Yq 7rfc c ln{6/a} 



b b 



expji/cr} 
cos <p dip dp dp , (9) 



where 



r = (p 2 + p' 2 - 2pp' coinp) 1/2 . 

i 2 2 

k = LU p,Q£o£ 



(10) 
(11) 



and k c is the wave vector of the coaxial line with the 
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characteristic admittance^ 



Y = 1/Zo = 2tt 



£q£ c 



ln{6/a} 



(12) 



The validity of the variational approximation ([9]) has 
been proven in Ref. Hi] in the parameter range < ka < 2 
and 1.57 < b/a < 4 by comparison with experimental 
results. To obtain low-temperature microwave data on 
Si:P with widely varying phosphorus concentration n, we 
employ a coaxial probe of dimensions 2b — 1.75 mm, 
2a = 0.6 mm. The maximum frequency range spans 
from 45 MHz to 40 GHz (limited by source and test set 
of the network analyzer HP 8510) and the relative dielec- 
tric constants £\ reach from 18 till 50. This corresponds 
to 0.001 < ka < 1.76 and b/a = 2.9 and lies within the 
tested parameter range. 



C. Solution of the inverse problem 

The inverse problem of extracting e from the measured 
impedance Z using Eq. j9| has been solved in the quasi- 
static approximation by Misraii. For low frequencies, 
the exponential function in Eq. ([9]) can be approximated 
by the first four terms of its series expansion: 
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(13) 



The second term of Eq. (|13[) vanishes upon integration, 
and the last one is readily integrated; the integrals cor- 
responding to the first and the third terms need to be 
numerically evaluated: 
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b b 



cos ip(p 2 + p' 2 — 2pp' cosp) 1 ^ 2 dpdpdp' . 



In our special case the relative contribution of the last 
term in Eq. (fT4)) to Y is below 1% up to 5 GHz and the 
formula thus reduces to a quadratic equation for e: 



1_ _ -i2cj£ a £ 
Z ~ ~ (ln{6/a}) 2 



U 2 p, I 3 £ Q £ 



(15) 



TABLE I: Geometrical integrals for the coaxial probe dimen- 
sions 2a — 0.6 mm and 2b — 1.75 mm 



Ii, mm 


7 3 , mm s 


r 4 

i4, mm 


7 5 , mm 5 


0.9084 


-0.2100 


-(tt/4) -0.4001 


-0.4047 



The values of the geometrical integrals for the special case 
of the coaxial probe with the inner and outer conductor 
diameters 2a = 0.6 mm and 2b = 1.75 mm are listed in 
Tab. HI It should be mentioned that the integrand of I\ 
diverges at p = p', ip — 0. That integral was numerically 
evaluated as the limit of a series of integrals n } which 
lower bounds ip n converge to p = 0. 



V. OPEN CALIBRATION STANDARD 

The frequency dependence of the open calibration 
standard with a known dielectric function £ can be ob- 
tained as follows. The expression ^ of the admittance 
Y as a function of e describes the open standard admit- 
tance correctly as long as the effect of the finite sample 
dimensions is negligible. Using a teflon block of the form 
shown in Fig. [3] and assuming its dielectric function to be 
£ = 2.03(1 + i 0.0002) in the GHz frequency ranged, the 
maximum electric field strength at the depth of 2 mm 
turns out to be far below 0.01 of its value at the sample 
surface for frequencies up to 10 GHz, so that the secon- 
dary reflections at the back side of the open standard can 
be neglected here. 

In order to obtain a closed expression Y opcn (uj) from 
the integral equation ([5]), the series expansion of the ex- 
ponential function can be used as in the previous section. 
In contrast to the inverse problem discussed in Sec. IIV1 
there is no need to spare at the accuracy truncating the 
series early here. The relative contribution of the sub- 
sequent term being far below 10 up to 10 GHz, the 
ultimate expression we use is: 



-^opo 



where 
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(16) 



h = -^(b 2 -a 2 ) 2 , 



b b 



h = J J J cosp(p 2 + p' - 2pp' cos ip) 3 ' 2 dp dp dp' . 

a a 

and k is defined in Eq. (jTTj) . The values of the geomet- 
rical integrals for the 2.4 mm coaxial probe are listed in 
Table H 

The frequency dependent reflection coefficient S^™ of 
the open calibration standard follows using Eq. ([2|): 
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The effect of the frequency dependence of Sif en on the 
conductivity <j\ and permittivity £% spectra compared 
to the dc assumption S°f cn — 1 is demonstrated in 
Figs. [5] and [5] on the example of two Si:P-samples with 
donor concentration n/n c of 0.56 and 0.9 relative to 
the concentration value at the metal-insulator transition, 
n c = 3.5 x 10 18 cm~ 3 . For samples with larger dielectric 
constant E\ and higher losses &\, as n/n c rises, the in- 
fluence of this correction slightly decreases. This is also 
what one would expect, when the electric properties of 
the material under investigation approach the metallic 
characteristics. 



VI. APPLICATION TO THE HOPPING 
TRANSPORT IN SI:P 

A. Dynamical conductivity 

With the method described in the previous parts, 
Sees. ITVl andlVl we study the frequency-dependent hop- 
ping transport in Si:P in order to explore the influence of 
electronic correlations 27 . At concentrations of phospho- 
rus in silicon below the critical value of n c = 3.5 x 10 18 
cm~ 3 , the donor electron states are strongly localized due 
to disorder in Anderson sens o 28 ' 29 . Since some degree of 
compensation by impurities of the opposite type is in- 
evitable, charge transport at low excitation energies is 
by variable-range hopping between the donor sites, ran- 
domly distributed in spac o 29 ' 30 . Thus, theoretical mo- 
dels for a disordered system with electron-electron inter- 
action are appropriate to interpret the electric conduc- 
tivity spectral. The main issue we address is that of 
power laws of the frequency-dependent conductivity at 
zero temperature: 

fii(w)— w Q , a = ai+ia 2 - (18) 

From the theory of resonant photon absorption by 
pairs of states, one of which is occupied by an elec- 
tron and the other one is empty, distinct limiting results 
are known for the conductivity power law depending on 
which of the relevant energy scales of the problem domi- 
nates over the others: 

Taking into account the Coulomb repulsion C/(r w ) if 
both states in a pair would be occupied by an electron 
(r u is the most probable hopping distance), Shklovskii 
and Efros derived ai(u) to be a sub-linear function of 
frequency, as long as the Coulomb interaction term do- 
minates over the photon energy 3 ^. At higher frequencies, 
in the opposite limit, the sub-quadratic behavior known 
from Mott^ - for non-interacting electrons is recovered. In 
addition, it is known that due to electronic correlations 
an area of reduced density of states is formed around 
the Fermi level, the so-called Coulomb gap A. For the 
conductivity of interacting electrons where the Coulomb 
term U(rJ) dominates over the photon energy but falls 
inside the Coulomb gap A, the reduction of the density 
of states leads to a stronger, slightly super-linear power 
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FIG. 8: (color online) Typical spectra of the measured real 
part of the conductivity for Si:P with relative donor concen- 
trations n/n c of 0.65, 0.83 and 0.90. The dashed lines are the 
fits by a two-parameter power law function. 



law 3 -^. In order to gain some insight into the effects 
of electronic correlations, it is required to extract tho- 
roughly the frequency-dependent conductivity and the 
related power-laws over a wide spectral range for a vari- 
ety of doping concentrations. 

In Fig. [S] the measured real part of the frequency- 
dependent conductivity is plotted on a log-log scale to 
identify the power law. The fits by a two-parameter func- 
tion <7i(w) = const. • uj a are shown by the dashed lines. 
The frequency dependence of the conductivity clearly fol- 
lows a super-linear power law in the whole doping range, 
where the exponent decreases with doping^ 7 .. From this 
we infer, that hopping transport takes place deep inside 
the Coulomb gap A. A super- linear conductivity power 
law was previously observed in Si:As and Si:P by Castner 
and collaborators 3 ^ 3 , using resonator techniques at cer- 
tain frequency within the range of the present work. Our 
results are in accord with the measurements on similar 
samples at higher frequencies (30 GHz to 3 THz) using 
optical technique o 34 ! 35 . In contrast, a sub- linear fre- 
quency dependence in the zero-phonon regime has been 
reported by Lee and Stutzmanr* 3 ^ based on experiments 
on Si:B in the microwave range and by Helgren et ali^ 7 . 
using quasi-optical experiments. 

As demonstrated in Fig. O the surface impedance ap- 
proach, mentioned in Sec. Ill Al yields a wrong (i.e. too 
strong) frequency dependence of the conductivity o\ (oj) 
for insulators. The evaluation of the impedance spectrum 
of a Si:P-sample is shown using the surface impedance 
formula ([4J in comparison to the solution of the equa- 
tion (|15p . The latter method yields a conductivity power 
law of approximately one, as expected from the theory 
outlined above. 
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FIG. 9: (color online) Frequency-dependent conductivity of 
a thick Si:P-sample with relative donor concentration of 65% 
measured at 1.6 K. In the first case (red line) the evaluation of 
the impedance data is done by the surface impedance concept 
(Sec. Ill A[) . Using the correct field distribution (Sec. IIV|I . 
a significantly different frequency dependence (blue line) is 
obtained from the same raw data. 



B. Dielectric function 

It is an additional advantage of a phase sensitive mea- 
surement to gain the dielectric function e\ from the ima- 
ginary part of the complex conductivity 4 : 



£i + ie 2 



1 



EqUJ 



(19) 



We denote by e the full complex dielectric function of 
Si:P, relative to the free space permittivity Eq. As the 
metal-insulator transition is approached upon doping n, 
the localization radius diverges^. As a consequence, 
the electronic contribution to the dielectric function is 
also expected to diverge following a power law when the 
metal-insulator transition is approached 38 : 



£l — £Si 



\l-n/n c \-< 



(20) 



where esi — H-7 is the dielectric constant of the host 
material Si. 

From our experiments^ we find that the dielectric 
function e\ is independent of frequency in the range from 
50 MHz to 10 GHz, taking the measurement uncertainty 
into account (Fig. [S]) . A fit with the function (f2"0"l) results 
in an exponent Q =0.71, as shown in Fig. [TOl 

In the framework of the effective medium approxima- 
tion C, — 1 is expected 38 . From the quasi-optical experi- 
ments on Si:P different results are reported. Helgren et 
ali^i observe a similar dependence of the values of the 
dielectric constant on the donor concentration (though 
uniformly shifted to lower values by 8). Hering et al. 35 
have observed values of £i as we have, but with a much 
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FIG. 10: (color online) Doping dependence of the low- 
temperature dielectric constant of Si:P. The solid line is the 
fit by the power law function (|20[) for the electronic dielectric 
function Ei — esi. 



stronger donor concentration dependence of the dielectric 
constant, resulting in a much higher exponent £ = 1.68. 
It is obvious that this enormous discrepancy calls for fur- 
ther experiments which are more accurate as far as this 
analysis is concerned. 



VII. CONCLUSIONS 

We have thoroughly analyzed the problem of extrac- 
ting the electrical conductivity and permittivity from the 
complex impedance measured in the microwave range 
using a network analyzer. While for thin metallic films 
and bulk metals simple relations are readily available, 
special care has to be taken in the case of semiconduc- 
tors and insulators where the electric field penetrates and 
decays over an appreciable distance. Already a static 
model with an approximate field configuration leads to 
reasonable results. Eventually we present a rigorous solu- 
tion of the problem with basically no restricting assump- 
tions. The advanced analysis is applied to the broadband 
impedance measurements of Si:P with different doping 
concentrations and spanning a wide range of frequency. 
The findings can now be compared to the theory and 
yield important insight into the effects of electronic cor- 
relations on the hopping transport at low temperature. 
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